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ABSTRACT 

We present detailed semi-analytical models for the formation of disk galaxies both in a Universe 
dominated by dark matter (DM), and in one for which the force law is given by modified Newtonian 
dynamics (MOND). We tune the models to fit the observed near-infrared Tully-Fisher (TF) relation, and 
compare numerous predictions of the resulting models with observations. The DM and MOND models 
are almost indistinguishable. They both yield gas mass fractions and dynamical mass-to-light ratios 
which are in good agreement with observations. Both models reproduce the narrow relation between 
global mass-to-light ratio and central surface brightness, and reveal a characteristic acceleration, contrary 
to claims that these relations are not predicted by DM models. Both models require SN feedback in order 
to reproduce the lack of high surface brightness dwarf galaxies. However, the introduction of feedback to 
the MOND models steepens the TF relation and increases the scatter, making MOND only marginally 
consistent with observations. The most serious problem for the DM models is their prediction of steep 
central rotation curves. However, the DM rotation curves are only slightly steeper than those of MOND, 
and are only marginally inconsistent with the poor resolution data on LSB galaxies. 

Subject headings: galaxies: formation — galaxies: fundamental parameters — galaxies: spiral — 
galaxies: kinematics and dynamics — galaxies: structure — dark matter. 



1. INTRODUCTION 

The most successful cosmological model to date con- 
sists of a universe dominated by cold dark matter (CDM), 
possibly with a non-zero cosmological constant. However, 
this paradigm faces several serious shortcomings. First is 
the high central densities predicted for virialized dark ha- 
los. Navarro, Frenk & White (1996, 1997) have claimed 
the existence of a universal density profile for virialized 
halos, with pdm °c r _1 at small radii (hereafter referred 
to as the NFW profile). More recently, higher resolu- 
tion simulations have argued that CDM halos should have 
even steeper cusps (Fukushige & Makino 1997; Moore et 
al. 1998; Moore et al. 1999b; but see Kravtsov et al. 1998). 
Several studies have argued that these steep cusps are in- 
consistent with the observed rotation curves of dark mat- 
ter dominated systems, such as dwarfs and low surface 
brightness galaxies (Flores & Primack 1994; Moore 1994; 
Burkert 1995; Navarro 1998; McGaugh & de Blok 1998a; 
Stil 1999). Secondly, Klypin et al. (1999) and Moore et 
al. (1999a) have indicated an additional problem for the 
CDM paradigm: simulations predict that galaxies should 
have many more satellite systems than observed. Al- 
though the reason for this large discrepancy may be related 
to observational biases, feedback from supernovae, or an 
ionizing background, it might also indicate that the CDM 
power spectrum has too much power on small scales. Both 
these problems may (partially) be solved with a mixture 
of hot and cold dark matter. However, the dark matter 
scenario is certainly starting to loose its appealing char- 
acter, if we have to rely on a hybrid mix of baryons, cold 
dark matter, hot dark matter, and a non-zero cosmological 
constant in order to explain the observations. The prob- 
lems with fine-tuning become more and more severe, and 
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it might therefore be worthwhile to explore some alterna- 
tives to the dark matter theory. 

Although the flat rotation curves of galaxies are gener- 
ally interpreted as evidence for the existence of dark mat- 
ter, they could also result from our theory of gravity (or 
inertia) being in error, such that the dynamical masses 
inferred from the usual Newtonian equations are overesti- 
mated. Several groups have tried to eliminate the need for 
dark matter by modifying the theory of gravity (e.g., Mil- 
grom 1983a; Sanders 1986; Mannheim & Kazanas 1989; 
Liboff 1992; Moffat & Sokolov 1996). So far, the modi- 
fied Newtonian dynamics (MOND) proposed by Milgrom 
(1983a) has been remarkably successful, and in this paper 
we compare several predictions from MOND with those 
based on dark matter (hereafter DM). 

Disk galaxies are, from a dynamical point of view, fairly 
simple systems in which the discrepancies between the lu- 
minous and the inferred dynamical masses (the "mass dis- 
crepancies") are least ambiguous. In addition, they have 
been extensively studied, with data available for systems 
that span several orders of magnitude in both mass and 
surface brightness. They are thus the ideal test-beds to 
compare predictions based on MOND and DM with ob- 
servations in the hope to discriminate between the two 
alternatives. 

For both MOND and DM to be considered viable theo- 
ries, they should be able to explain the following properties 
of disk galaxies: 

1. The slope, scatter, and zero-point of the Tully-Fisher 
relation (Tully & Fisher 1977). In particular, the 
fact that both high and low surface brightness disks 
(HSBs and LSBs, respectively) follow the same re- 
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lation without any significant offset (Sprayberry et 
al. 1995; Zwaan et al. 1995; Hoffman et al. 1996; 
Tully & Verheijen 1997). 

2. The shape of the observed rotation curves. In par- 
ticular, the systematic variation of these curves with 
surface brightness and luminosity (Rubin, Thonnard 
& Ford 1980; Burstein et al. 1982; Casertano & van 
Gorkum 1991; Persic, Salucci & Stel 1996). 

3. The increase of mass discrepancies with decreasing 
luminosity and surface brightness (e.g., Salucci & 
Frenk 1989; Persic & Salucci 1988, 1990; Kormendy 
1990; Broeils 1992; Forbes 1992; Persic, Salucci & 
Stel 1996; McGaugh & de Blok 1998a). This is di- 
rectly linked to item 2. 

4. The increase of gas mass fractions with decreas- 
ing luminosity and surface brightness (Gavazzi 1993; 
Gavazzi, Pierini & Boselli 1996; McGaugh & de Blok 
1997). 

5. The sharp cut-off in the distribution of central sur- 
face brightnesses at the bright end, close to the value 
of the Freeman (1970) law (Allen & Shu 1979; Mc- 
Gaugh, Bothun & Schombert 1995; de Jong 1996a; 
McGaugh & de Blok 1998a). 

6. The presence of a characteristic acceleration in the 
rotation curves of disk galaxies (McGaugh 1998). 

In a series of papers, McGaugh & de Blok (1998a, b; 
hereafter MB98a and MB98b) have used several of these 
observational facts to argue against the DM hypothesis 
and in favor of MOND. In particular, they argue that 
the DM hypothesis faces severe fine-tuning problems. In 
this paper we revisit these issues using more sophisticated 
models for the formation of disk galaxies under both the 
DM and MOND hypotheses. We show that, contrary to 
the claims made by MB98a and MB98b, the DM model 
can explain virtually all aforementioned observations with 
a minimal amount of fine-tuning. Furthermore, we show 
that MOND needs a similar amount of fine-tuning, there- 
with loosing its apparent advantage over the DM scenario. 

This paper is organized as follows. In ^ we discuss 
the models, both for the DM and the MOND models. In 
§[| we use the near-infrared Tully-Fisher relation to con- 
strain the free parameters of the models. In §|] we com- 
pare the models to data on gas mass fractions, dynamical 
mass-to-light ratios, and characteristic accelerations. §|f] 
compares the detailed rotation curve shapes from the DM 
and MOND models. We summarize our conclusions in 
§||. Throughout this paper we adopt a Hubble constant of 
H = 70h 70 kms _1 Mpc _1 in all our models. Any quan- 
tity that depends on the distance scale is written in terms 
of h m . 

2. THE FORMATION OF DISK GALAXIES 

2.1. The dark matter hypothesis 

In van den Bosch (1999; hereafter Paper I) we presented 
detailed models for the formation of disk galaxies in the 
context of a universe dominated by CDM. Here we briefly 
summarize the ingredients of the models, and refer the 
reader to Paper I for details. 



In our models, disks form by the settling of baryonic 
matter in virialized dark halos, which are described by the 
NFW density profile. We assume that while the baryons 
radiate their binding energy, they conserve their specific 
angular momentum, thus settling into a disk with a scale 
length that is proportional to the angular momentum and 
size of the dark halo (see also Fall & Efstathiou 1980; 
Dalcanton Spergel & Summers 1997; Mo, Mao & White 
1998; van den Bosch 1998). The angular momentum, J, 
is thought to originate from cosmological torques and can 
be characterized by a dimensionless spin parameter, A, the 
distribution of which is well constrained from both analyt- 
ical and numerical studies (e.g., Barnes & Efstathiou 1987; 
Warren et al. 1992). For each separate model galaxy, we 
randomly draw a value for A from this distribution. The 
disk is assumed to be a thin exponential with surface den- 
sity S(i?) = So exp(— R/Rd). Adiabatic contraction is 
taken into account to describe how the density profile of 
the dark matter changes due to the cooling of the bary- 
onic matter (Blumenthal et al. 1986; Flores et al. 1993). A 
global stability criterion is used to investigate if the disks 
that form are stable. If unstable, part of the disk mate- 
rial is transformed into a bulge component, until the disk 
reaches stability (van den Bosch 1998). 

Once the density distribution of the baryonic material is 
known, we compute the fraction of baryons converted into 
stars. Kennicutt (1989) has shown that the star formation 
rate (SFR) in disk galaxies follows a Schmidt law (Schmidt 
1959), but is abruptly suppressed below a given threshold 
density. The critical density is given by Toomre's (1964) 
stability criterion as: 
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Here Q is a dimensionless constant near unity, cr gas is the 
velocity dispersion of the gas, and k is the epicycle fre- 
quency, which can be derived from the rotation curve. 

Setting dgas = 6 kms -1 and Q = 1.5 these star for- 
mation thresholds were found to coincide with the radii 
where E gas = S cr it. Following an earlier suggestion by 
Quirk (1972), we therefore consider the fraction of disk 
mass with surface densities in excess of S cr i t to be eligible 
for star formation. If one further considers the empirical 
Schmidt law determined by Kennicutt (1998) for a large 
sample of disk galaxies, it is found that over the typical 
lifetime of a galaxy virtually all the mass which is eligi- 
ble for star formation does in fact turn into stars. This 
is consistent with the lack of major star formation in the 
inter-arm regions in spirals and yields gas mass fractions 
that are in good agreement with observations (Paper I). 
In what follows we therefore assume that the entire bulge 
mass and the entire disk mass with S(i?) > £ C rit are trans- 
formed into stars, giving a total stellar mass of M». 

The total if -band luminosity of the disk-bulge system is 
then calculated using Lk — M*/T* K , where is the stel- 
lar mass-to-light ratio. Since depends only weakly on 
the age and metallicity of a stellar population (see discus- 
sion in Paper I), we assume a constant stellar mass-to- light 
ratio of = 0.4 h 70 M / L Q . 

With the mass in stars known, the total amount of en- 
ergy produced by SN can be calculated, once we make 
assumptions about the energy per SN (-Esn) and the num- 
ber of SN per solar mass of stars formed (?7sn)- We then 
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calculate the total baryonic mass that is prevented from 
becoming part of the disk/bulge system due to SN feed- 
back as 

M hot = 3.22M*eg N x 

( vs-n \ ( e sn \ i v 200 y- 2 ( 9 \ 

\^O.OO4M - i y \^10 51 cr g y V250kms- 1 / ' \") 

Here V200 corresponds to the circular velocity of the halo 
at the radius r2oo, inside of which the average density 
equals 200 times the critical density for closure, and £g N 
and v are two free parameters that describe the efficiency 
with which the SN energy is used to prevent gas from 
settling in the disk/bulge system (see Paper I). Through- 
out we adopt a value of ??sn = 0.004 M Q _1 which cor- 
responds to a Scalo (1986) initial mass function. An it- 
erative procedure is used to compute the masses in the 
various phases in a self-consistent manner, i.e., we ensure 
that M 200 = A/dm + M* + M coW + M hot , with Af 20 o the 
total mass within r2oo (before SN blow-out), and M co \d 
the disk mass that is not turned into stars. 

2.2. The modified Newtonian dynamics hypothesis 

We now construct similar models but under the hypoth- 
esis of MOND, which assumes that the force law changes 
from conventional Newtonian form when the acceleration 
of a test particle is smaller than a characteristic accelera- 
tion do (which is a universal constant). The true gravita- 
tional acceleration a is related to the Newtonian accelera- 
tion aN as 

^(a/a )a = a N (3) 

(a = |a|). The interpolation function fi{x) is not spe- 
cific, but is required to have the asymptotic behavior 
fi{x » 1) — ► 1 (the Newtonian regime) and [i(x <C 1) — * x 
(the MOND regime). With the rotation law given as usual, 
U c 2 (r)/r — a, the circular velocity of a finite bounded mass 
M in the limit a — ► becomes 



Voo = (GMa ) 1/4 , 



(4) 



with G the gravitational constant. MOND thus predicts 
asymptotically flat rotation curves. 

First we need to set ao and the interpolation func- 
tion /i(x). We follow previous studies of MOND, and set 
a = 1.12 x 10~ 10 h 70 ms -2 and adopt 

fi(x) = 



(5) 



For these choices, MOND yields good fits to the observed 
rotation curves of disk galaxies (Kent 1987; Milgrom 1988; 
Begeman, Broeils & Sanders 1991; Sanders 1996; de Blok 
& McGaugh 1988). With fi{x) given by equation (§), the 
circular velocity as function of radius becomes 



V c (r) = V N (r) 




1/4 



(6) 



where Vjv(r) = yr d<I>/dr is the circular velocity derived 
for pure Newtonian dynamics. 

Unfortunately, there is currently no MOND equivalent 
to the theory of structure formation in a DM universe. 



In particular, we don't know the formation epochs or dis- 
tribution of angular momenta of galaxies that form in a 
MOND universe. We therefore follow an empirical pro- 
cedure in which we assume that galaxy formation under 
MOND produces disk galaxies with the masses and di- 
mensions as observed. As in the DM scenario, disks are 
assumed to be thin exponentials. We randomly draw a to- 
tal (baryonic) mass from the interval lO^^M© < M to t < 
5.7 x lO 11 /).^ 1 M . We initially set the disk mass, M<z, 
equal to M to t, and randomly draw a central disk sur- 
face density from the interval 45/i7o M pc~ 2 < So < 
2400 hro M Q pc~ 2 . The scale length of the disk is then 

i? d = (M d /2 7rI]o) 1/2 - 

The first step is to investigate whether the disk is sta- 
ble against global instabilities. To that extent we use the 
same stability criterion as for the DM scenario, i.e., a disk 
is considered stable if 



°' 3 V c (3R d ) < a ^ 



(7) 



where V c (3Rd) is the circular velocity of the galaxy at three 
scale lengths (cf. Christodolou, Shloshman & Tohlinc 
1995; Efstathiou, Lake & Negroponte 1982; van den Bosch 
1998). In the case of Newtonian dynamics with no DM ha- 
los, a — 0.5. However, in the case of MOND, a < 0.5, at 
least if the acceleration at 3 Rd is sufficiently small; i.e., the 
MOND force law helps to stabilize disks (Milgrom 1989). 
Unlike the DM scenario, for which numerical simulations 
have shown that a cr it — 0.35 (Christodolou et al. 1995), no 
equivalent study exists based upon MOND. We therefore 
consider a cr it a free parameter. 

Once the density distribution of the gas is known, we 
calculate the mass in stars, in the same fashion as for the 
dark matter models discussed above. Note, however, that 
in the calculation of the epicycle frequency, k, we now use 
the circular velocity based on the modified force law. The 
i^-band luminosity is again calculated from M* adopting a 
constant mass-to-light ratio , which we consider a free 
parameter. 

Finally, we include SN feedback, using the same ap- 
proach as above for DM. The hot gas mass, Mhot is cal- 
culated using equation (||), but with V200 replaced by Voo 
(equation j4|). We use an iterative technique to ensure 
self-consistency, i.e., M tot = M* + M cold + Af ho t- 

3. CONSTRAINTS FROM THE TULLY-FISHER RELATION 

The slope, scatter and zero-point of the Tully-Fishcr 
(hereafter TF) relation contain information about the 
structure and formation of disk galaxies. Here we use the 
observed characteristics of this fundamental scaling rela- 
tion to set the parameters of our disk formation models 
described in the previous section. 

3.1. Observations and simple dynamical predictions 

The slope and scatter of the empirical TF relation de- 
pend strongly on the luminosity and velocity measures 
used, with the relation becoming shallower going towards 
bluer passbands (Aaronson, Huchra & Mould 1979; Vis- 
vanathan 1981; Tully, Mould & Aaronson 1982; Wyse 
1982; Pierce & Tully 1988; Gavazzi 1993; Verheijen 1997). 
Therefore, when trying to infer constraints on the struc- 
ture and formation of disk galaxies from the TF relation, it 
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is absolutely essential that one extracts the same luminos- 
ity and velocity measures from the models (see discussions 
in Courteau 1997, Verheijen 1997, and Paper I). 

The empirical TF relation which most directly reflects 
the mass in stars and the total dynamical mass of the halo 
is the i-T-band TF relation of Verheijen (1997), which uses 
the flat part of the rotation curve as a characteristic ve- 
locity: 



L K = 2.96 x 10 11 L c . 



V 6a 



250 km s" 



4.2 



(8) 



The slope of this TF relation is b = —10.5 ± 0.5, and the 
observed scatter is <jm — 0.29 mag. This empirical TF 
relation is shown in the upper left panel of Figure [j]. Open 
circles correspond to the galaxies in Verheijen's sample of 
Ursa Major spirals, and the thick solid line is the best- 
fitting TF relation (equation ||]). 

Within the DM framework, simple dynamics predict a 
TF relation of the form 



L 



K 



e gf fba 



10GH T K VVa. 



flat' 



(9) 



(Dalcanton et al. 1997; White 1997; Mo et al. 1998; Syer, 
Mao & Mo 1999; van den Bosch 1998; Paper I). Here e gf 
is the galaxy formation efficiency, which describes what 
fraction of the baryonic mass inside the DM halo ulti- 
mately ends up in the disk/bulge system, /bar is the bary- 
onic mass fraction of the Universe (/bar = ^bar/^o), and 
Tk = Mgai/Lx = T^Mgai/M* is the if-band mass-to- 
light ratio of the galaxy (disk plus bulge). 

From equation (jij) it is immediately apparent that in 
the case of MOND, the TF relation becomes 



J K 



cioGYk \Vfi 



V, 



flat J 



(10) 



(see e.g., Milgrom 1983b; MB98b). 

Thus whereas DM predicts a TF relation with a slope 
of b = —7.5, MOND predicts b — —10, in agreement with 
the empirical relation of equation (||). In the case of DM, 
extra physics are required to steepen the predicted TF 
relation to the observed slope. MB98a have presented a 
detailed discussion in which they argue that this results in 
a serious fine-tuning problem for DM. On the other hand, 
the almost perfect agreement between the observed TF 
slope and the pure-gravitational value for MOND (equa- 
tion |To| ) presents an anti-fine-tuning problem. There is 
very little room for any other galaxy characteristic to vary 
systematically with mass, without changing the TF slope. 

The above suggests that simple dynamical arguments 
can be used to predict the existence of an underlying TF 
relation for both DM and MOND. However, to accurately 
predict the slope, zero-point, and scatter of the TF rela- 
tionship, it is necessary to incorporate the more detailed 
physics discussed in §0. We do so in the following two sec- 
tions for the D M a nd MOND hypotheses, respectively. As 
mentioned in §3.1 it is essential to extract the same veloc- 



ity measure from the models as the one used in the data, 
and we therefore extract Vg a t from our model galaxies (see 
paper I for details). 



3.2. The Tully-Fisher relation and dark matter 

In Paper I we showed that taking the stability-derived 
star formation threshold densities into account (with Q = 
1.5) steepens the TF relation from Lk oc Vg at (as pre- 
dicted by pure dynamics) to Lk oc V^. We argued that 
to further steepen the TF relation to its observed slope, 
feedback from SN is required. Fine-tuning the two param- 
eters that describe the efficiency with which SN energy is 
used to prevent gas from settling in the disk/bulge system 
(eg N and is), we were able to obtain a TF relation with 
a slope, zero-point, and an amount of scatter that are all 
in excellent agreement with the empirical relation of equa- 
tion (||) . This model is referred to as model L5 in Paper I, 
and we adopt the same convention here. This model is 
based on a cosmology with f^o = 0.3, f^A = 0.7, hfo = 1 
and erg = 1.0 (see paper I), and its parameters are listed 
in Table [l| Model L5 predicts an intrinsic TF scatter of 
o~m ~ 0.2 mag, with no offset in normalization between 
HSB and LSB galaxies, as observed. The TF relation of 
model L5 is shown in the upper right panel of Figure [l| 

3.3. The Tully-Fisher relation and MOND 

Although simple dynamics predict a TF relation for 
MOND with a slope b = —10, consistent with observa- 
tions, this ignores the fact that Tk varies with mass; gas 
mass frac tions are known to increase with decreasing mass 
(see §p~2|), and taking this into account steepens the TF 
relation. It therefore remains to be seen whether MOND 
does not actually predict TF relations that are too steep. 



We can investigate this using the models described in §2.2 
As with the DM models, we set Q = 1.5. The stel- 
lar mass-to- light ratio in the if-band, T K , is set by fit- 
ting the normalization of the TF relation, yielding T K — 
0.53 /170 Mq/L©. This value is consistent with estimates 
based on stellar population models, which in itself can be 
considered a success of MOND. 

To set the free parameter a cr it, with the help of equa- 
tion (^) the stability criterion (equation can be written 



S < 1.91 x 10 3 (- 



«o 



-1 -1 



-1/2 



M pc- 2 .(11) 



The stability parameter a cr it thus directly translates into a 
maximum surface density of disks. From the fact that the 
brightest observed disks have central surface brightnesses 
of /xo,k — 16-5 mag arcsec -2 (e.g., de Jong 1996b; Tully et 
al. 1996) we infer a crit > 0.45 (for = 0.53 h 70 M /L Q ). 
If we were to take a cr it = 0.35 (as for the DM models), 
no disks with central surface brightnesses brighter than 
A*o, k = 18.1 mag arcsec -2 would be allowed, clearly in- 
consistent with observations. Since in general not all the 
gas mass is transformed into stars, the limit of a cr i t = 0.45 
is conservative. We find that with Q = 1.5 and T K = 
0.53 hro M Q /L a value of a cr it = 0.48 yields /io,ir <; 16.5 
mag arcsec -2 , and we adopt this value throughout. None 
of the results presented here, however, depend on this 
choice: a cr it merely sets the maximum central surface den- 
sity for disks with a given bulge-to-disk ratio, and is set 
such that we obtain disks with the full range of observed 
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values of [Iq^k- In mimicking the sample selection of Ver- 
heijen (1997) we only consider model galaxies with a bulge- 
to-disk ratio B/D < 0.2 (see also Paper I). 

We compare the DM model (L5) with two different 
MOND models, Ml and M2. The former neglects SN feed- 
back, whereas the latter considers blow-out by SN with 
eg N = 0.05 and v = -3.0 (see Table |). The motiva- 
tion for t hese particular feedback parameters will be dis- 
cussed in §[D^. The TF relations for models Ml and M2 are 
shown in the lower panels of Figure [j]. The star formation 
recipe with Q = 1.5 steepens the MOND TF relation from 



L K oc U fl 4 at (equation (10|) to L K oc U fl 4 at 2 (model Ml) 



and brings it in excellent agreement with the empirical 
TF relation of Verheijen. The predicted, intrinsic amount 
of scatter is only <tm — 0.21 mag, once again in excel- 
lent agreement with observations. The introduction of SN 
feedback steepens the TF relation further, to b = —11.9 
for model M2, and increases the scatter to <jm — 0-33 mag. 
Model M2 thus yields a TF relation that is too steep. In 
addition, the amount of scatter is only marginally consis- 
tent with the data, i.e., Verheijen (1997) derived <jm < 0.3 
mag at 95 percent confidence level. 

4. COMPARISON OF MODELS WITH DATA 

Having tuned the parameters by fitting the near-infrared 
TF relation, we now compare the resulting models to other 
independent observations, to test, amongst others, the va- 
lidity of the star formation and feedback recipes used. As 
shown in paper I, within the CDM context one can con- 
struct many physically different models that all fit the 
empirical TF relation. For these models additional con- 
straints, such as the ones described below, are required to 
descriminate between the models. 

4.1. The data 

Our calculations yield magnitudes and surface bright- 
nesses in the i^-band and rotation velocities for a collec- 
tion of model galaxies. Ideally, we would therefore com- 
pare our models to data that consists of -ftT-band photome- 
try combined with full HI rotation curves. Unfortunately, 
such data is currently scarce, and we therefore use J3-band 
photometry, for which ample data is available. Although 
a constant stellar mass-to-light ratio is a reasonable as- 
sumption for the X-band, T B is expected to reveal a sys- 
tematic trend with mass as well as a significant amount of 
scatter owing to variations in dust content, and the ages 
and metallicitics of the stellar populations. We convert the 
X-band luminosities and surface brightnesses of our model 
galaxies to the -B-band adopting an empirically determined 
B — K color magnitude relation (see Appendix A). This, at 
least, takes account of the systematic trend of B — K with 
Mk- Unfortunately, the scatter is relatively large, and we 
predict errors in Mb and Hq^b, owing to the color conver- 
sion, of 0.37 and 0.47 mag, respectively. A more detailed 
comparison of our models against data has to await ac- 
curate near-infrared photometry of a large sample of disk 
galaxies for which full HI rotation curves are available. 

We compiled the following data for comparison with our 
models: 

• Mk, (J-o.k, Rd, and Vflat for a sample of spirals in 
the Ursa Major cluster from the thesis of Verhei- 
jen (1997). We restrict ourselves to the twenty-two 



galaxies of the "unperturbed sample" on which the 
TF relation of equation (||) is based. As for our mod- 
els, we transform Mk and (Iq.k to the -B-band, using 
the B — K color magnitude relation presented in Ap- 
pendix A. Data from this sample are plotted as open 
circles. 

• Mb, Mo,-Bi Rd, and Ufl a t for a sample of thirty-two 
galaxies compiled by MB98a. These data are plotted 
as solid circles. 

• Mb, Ho,b, Rd, and Vm ax for ten dwarf galaxies from 
the sample of van Zee et al. (1997). The velocity 
measure U max is virtually identical to the velocity at 
the last measured point of the HI rotation curves, 
and as such compatible with the velocity measure 
determined in our models (see paper I). These data 
are plotted as solid squares. 

• Mb, Mo.-Bj and Mm/ Lb for a sample of one hundred 
disk galaxies of type Sb or later, compiled by Mc- 
Gaugh & de Blok (1997). These data are plotted as 
open squares. 

All data have been converted to hjo = 1, using the dis- 
tance scales quoted by the various authors. Note that this 
combined set of data is by no means complete nor homo- 
geneous. In compiling the data we have aimed for a max- 
imum coverage in both luminosity and surface brightness 
of disk galaxies with measured HI rotation curves. There- 
fore, it is important to keep in mind that upon comparing 
models to data, one should only focus on global properties, 
such as the areas of parameter space occupied by the data 
and the models. Our models are not constructed to cor- 
rectly predict the local density distribution of data points 
in parameter space (i.e., the masses of our model galaxies 
are sampled uniformly rather than from a Press-Schechter 
formalism or a luminosity function). In addition, the in- 
completeness and possible observational biases in the data 
samples, combined with the errors associated with the con- 
version of X-band models to the -B-band, prevents us from 
testing our models in too much detail. 

4.2. Gas mass fractions 

One of the empirical characteristics of disk galaxies is 
the increase of their gas mass fractions with decreasing lu- 
ll). This is 
, where we 



minosity and surface brightness (i.e., item 4 in 
shown in the upper left panels of Figures ^| and 
plot the logarithm of the HI mass-to-light ratio, Mm/ Lb, 
as function of Mb and Ho,b, respectively. The remaining 
panels plot the results for each of the three models. The 
solid lines in Figures [2] and || have no physical meaning, 
but are plotted to facilitate a comparison of the models 
with the data. 

The increase of Mhi/Lb with decreasing magnitude 
and surface brightness is remarkably well reproduced by 
each of the three models. Note also that the range spanned 
by the model galaxies is in good agreement with the data. 
The HI mass-to-light ratios derived from the models scale 
with /170, whereas the measured ratios are independent of 
distance (the HI mass is computed directly from the ob- 
served flux). The main parameter in the models that sets 
Mhi/Lb is the star formation parameter Q. It is reassur- 
ing that for H = 70kms~ 1 Mpc~ x and for the value of Q 
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that was empirically determined by Kennicutt (1989), the 
models yield gas mass fractions that are in good agreement 
with observations. 

4.3. Dynamical mass-to-light ratios. 

Dynamical mass-to-light ratios are seen to vary system- 
atically with mass and surface brightness. Dwarfs and LSB 
galaxies have consistently higher mass-to-light ratios than 
bright HSB galaxies. These mass-to- light ratios are higher 
than can be explained by differences in stellar populations 
alone, and indicate an increasing "mass discrepancy" with 
decreasing luminosity and surface brightness (i.e., item 3 
in§§. _ 

In this paper we define a characteristic mass-to-light ra- 



tio^ 



T 



R dVj at /G 

L B 



(12) 



In Figure [| we plot log(To) as function of Mb for the 
data and our three models. The thin line is plotted to fa- 
cilitate a comparison, and is chosen to roughly outline the 
region that reveals a clear deficit of galaxies; low luminos- 
ity galaxies always have high mass-discrepancies, whereas 
bright disk galaxies can have values of To that span over 
an order of magnitude. This absence of faint galaxies with 
low mass discrepancies is well reproduced by model L5, 
and owes largely to the effect of SN feedback. The effi- 
ciency with which SNe can blow gas out of a dark halo 
is larger in less massive galaxies. This causes most of the 
baryonic mass to be expelled from low-mass halos, result- 
ing in higher values of To for fainter galaxies. A possible 
shortcoming of model L5 seems to be the fact that it pre- 
dicts virtually no galaxies with T < 0.8, whereas four of 
the sixty-four galaxies in the combined data set have val- 
ues of To below this value. However, given the scatter in 
the B — K color magnitude relation, the inhomogeneous 
nature of our comparison sample, and the typical errors 
associated with the data, this apparent deficit is not sig- 
nificant. 

In the MOND model without feedback (model Ml), 
there is no reason for fainter galaxies to have systemat- 
ically higher values of To- Indeed model Ml does not 
reproduce the observed absence of dwarf galaxies with low 
mass discrepancies, and it does a poor job in reproducing 
the data. We therefore included feedback, and tuned the 
parameters £g N and v such that the observed 'region of 
avoidance' is reproduced, yielding model M2. Note that 
model M2 predicts, similarly to model L5, an absence of 
model galaxies with To <J 0.8. 

A similar result is shown in Figure ^, where we plot 
the .B-band central surface brightness versus Va a t- The 
data reveal a clear absence of HSB disk galaxies with low 
rotation velocities. In other words, disk galaxies with 
Vflat <J 50 kms -1 always have low surface brightness. 
Again, our DM model reproduces this behavior remark- 
ably well. Model Ml predicts no such deficit, and is in 
clear conflict with the data. Model M2, however, in which 
the two SN parameters have been tuned to yield the ob- 
served absence of dwarf galaxies with low mass discrepan- 
cies, is in excellent agreement with the data. 

Note that the observed regions of avoidance in both the 



Mb - To and the Vg a t - /io,B planes are not due to ob- 
servational biases: There is no reason why dwarf galaxies 
with high central surface brightnesses should be missed, 
whereas their LSB counterparts are not (except for a pos- 
sible problem due to star-galaxy separation for compact 
systems at large distances). This therefore seems to be 
a strong indication that SN feedback plays an important 
role in the formation of (disk) galaxies, at least for the less 
massive ones. 

MB98a have shown that To is strongly correlated with 
central surface brightness. This is confirmed from our com- 
bined data set (which includes the data of MB98a) plotted 
in the upper left panel of Figure ^. The thin solid line cor- 
responds to Tq oc 1/Eq, ana - is plotted with arbitrary nor- 
malization for comparison. Here Eg is the central surface 
brightness of the luminous, stellar disk in L /pc 2 . The 
data clearly indicates that Tg Eg ~ constant with only 
little scatter. 

MB98a and MB98b argue that this To-Eo "conspiracy" 
follows naturally from MOND, whereas it requires a prob- 
lematic fine-tuning in the case of DM. As shown by Zwaan 
et al. (1995), 

To ^o = constant is required to obtain a TF 
relation in which HSB and LSB galaxies follow the same 
relation with a slope of b = —10. MB98a discuss a num- 
ber of physical processes that may lead to the observed 
Tg -Eg "conspiracy" under the hypothesis of DM, none 
of which they render feasible. The upper right panel of 
Figure [|, however, shows that our DM model reproduces 
the observed "conspiracy" to a good degree of accuracy. 
Given that the same model yields a TF relation with a 
slope close to —10 and with little scatter, this is not too 
surprising, since the To~Eo conspiracy and the TF relation 
are closely connected (see discussions in Zwaan et al. 1995 
and MB98a). We thus disagree with MB98a, and conclude 
that the To~Eo relation is in fact fairly easily reproduced. 

The lower two panels of Figure ^ show that the two 
MOND models are also in good agreement with the data, 
something which is not too surprising given that pure 
MOND dynamics predict that T§ cx l/Y,* (MB98b). 



4.4. Characteristic 



derations 



Milgrom (1983b) defined the characteristic acceleration 
parameter 



V? 



flat 



ao Rd 



(13) 



which is proportional to the ratio of the characteristic ac- 
celeration Vfl^/Rd of a disk and the characteristic MOND 
acceleration ao. As shown by MB98b, for an exponential 
disk MOND predicts a tight correlation between £ and 

surface brightness of the form £ cx Eg. The upper left 
panel of Figure M plots log(£) as function of /io.b for the 

. ' 1/2 

data. The thin solid line corresponds to £ oc E , and 
is plotted with an arbitrary normalization. Clearly, the 
data reveals a narrow correlation between £ and central 
surface brightness, with a slope close to that predicted 
by MOND. MB98b considered this yet another success of 
modified Newtonian dynamics. 

However, as we show in Appendix B, the dark matter 

1/2 

hypothesis predicts the same relation, i.e., £ oc E , and 

and that disks are 



2 MB98a denned To as the global mass-to-light ratio inside four disk scale lengths. If we assume that Vc('iRd) — Vg 
pure exponentials, our definition of To is identical to that of MB98a, apart from a factor 0.23. 
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the observed relation between £ and (Iq.b thus does not 
provide a useful test to discriminate between the MOND 
and DM hypotheses. This is also evident from Figure 0, 
which shows that each of our three models yields a relation 
between £ and fJ-o,B that is in excellent agreement with 
the observations. It is important, though, to note that 
under the DM hypothesis, the amount of scatter can be 
significantly larger than observed. It is because of the SN 
feedback, the efficiency of which decreases with increasing 
halo mass, that the amount of scatter in model L5 is as 
small as it is (see Appendix B for a detailed discussion). 

Additional support for MOND was recently presented 
by McGaugh (1998), who showed that current observa- 
tions hint to the presence of a characteristic acceleration, 
the main assumption on which MOND is based. The up- 
per panels of Figure || plot data, kindly provided by Stacy 
McGaugh, of the cumulative mass discrepancy T against 
radius R (panels on the right), against orbital frequency 
w = V I0 t/R (middle panels), and against centripetal ac- 
celeration a = V^ ot /R (panels on the right). Here T is 
defined as 



T(i?) = 



RV T l t (R) 



G[M*(R)+M coU (R)]' 



(14) 



and corresponds to the ratio of total mass (based on the 
Newtonian equations) to luminous mass (stars plus cold 
gas) within radius R. Each data point represents one re- 
solved measurement in the rotation curve of a disk galaxy. 
The galaxies in the sample span a wide range in lumi- 
nosities and surface brightnesses and were taken from the 
compilation of Sanders (1996) and de Blok & McGaugh 
(1998). 

The plot of T versus R reveals that the mass discrep- 
ancy sets in at a wide range of different radii; individual 
rotation curves are easily discerned. This scatter is re- 
duced when plotting T versus the circular frequency, and 
becomes minimal when plotted versus the acceleration a. 
The data clearly reveal that in each disk galaxy, the mass 
discrepancy starts to become significant below an acceler- 
ation of ~ 10~ 10 ms~ 2 , as predicted by MOND. 

The lower two rows of panels of Figure || show the same 
plots for the two MOND models. Here we have randomly 
selected 40 galaxies from models Ml and M2, and for each 
model galaxy, we computed T at 15 radii, sampled uni- 
formly between R = and the radius at which the pro- 
jected HI column density is equal to 10 20 cm -2 . The two 
MOND models reproduce the data extremely well, which 
is to be expected, given the presence of the characteristic 
acceleration ao and the fact that MOND was specifically 
designed to fit rotation curves. More remarkably, model 
L5 also yields a characteristic acceleration (see panels in 
second row of Figure ||). Although the agreement with 
the data is not quite as good as in the case of models Ml 
and M2, the main characteristic of the data, namely the 
minimization of scatter when T is plotted versus a, is re- 
produced well. Most importantly, this is achieved without 
any fine-tuning of our model; after setting the feedback pa- 
rameters to fit the slope of the empirical TF relation, the 
same model reveals a characteristic acceleration. This is 
in clear contradiction with McGaugh (1998), who argued 
that this phenomenology can not be reproduced with DM 
models. Not only is this a remarkable success of our DM 



model, it is also an argument against an argument in favor 
of MOND. 

Model L5 predicts mass discrepancies which are slightly 
too large in the range 3x 10 9 < a < 10 10 ms~ 2 . The origin 
of this small deficit is discussed in the following §. 

4.5. Understanding the scaling relations 

Why can our DM models easily account for the To-Eo 
"conspiracy", whereas MB98a have argued that no feasi- 
ble solution exists? Combining equations |l^] and |l| , and 
taking into account that the total luminosity of an expo- 
nential disk scales as L cx Eo-R 2 ), one obtains To oc £ Eg -1 . 

1/2 

Since, as shown in Appendix B, we have £ cx E one thus 



automatically predicts that Tq oc E, 



-1/2. 



i.e., there is no 



conspiracy. However, the fact that the scatter around the 
To-Eo relation is so small is not trivial, as the scatter in 
the £-£o relation can be fairly large. It is only because 
of our particular feedback model, which we tuned to re- 
produce the empirical TF relation, that the scatter in the 
£-Eo relation, and therewith in the To~Eo relation, is so 
small (see Appendix B). We consider this another major 
succes for our feedback model. 

From the scaling relations listed above one derives that 
£ cx Tq 1 . We thus expect a narrow correlation between 
the characteristic accelerations and mass-to-light ratios of 
disk galaxies. Therefore, the fact that our DM model is 
succesful in reproducing the narrow correlation between 
local accelerations and mass-to-light ratios is not too sur- 
prising. 

5. ROTATION CURVE SHAPES 

The rotation curves (hereafter RCs) of disk galaxies are 
the most direct indicators of a discrepancy between the 
luminous mass and the inferred dynamical mass for disk 
galaxies. The shapes of rotation curves therefore provide 
a clean laboratory for comparing models of gravitational 
dynamics. Historically, astronomers have explained the 
detailed RCs of spiral galaxies by invoking the presence 
of non-luminous matter (see reviews in Sancisi & van Al- 
bada 1987, and Ashman 1992), while preserving the r 
force law of Newtonian gravity. More recently, some as- 
tronomers have used MOND to explain the observed kine- 
matics without requiring the presence of dark matter (e.g. 
Kent 1987; Lake 1989; Begeman et al. 1991; Sanders 1996; 
de Blok & McGaugh 1998; Sanders & Verheijen 1998). 
Both approaches have found nearly equal success over a 
wide range of luminosities and surface brightnesses, well 
beyond the luminous, high-surface brightness galaxies for 
which the theories of DM and MOND were developed. 

Recently, an important problem for the CDM picture 
has emerged from HI studies of dwarf galaxies. The slowly 
rising rotation curves observed in these systems suggest 
a (close to) constant density core in their dark matter 
halos, inconsistent with the steeply cusped density pro- 
files predicted for CDM (Flores & Primack 1994; Moore 
1994; Burkert 1995; Burkert & Silk 1997; Stil 1999, but see 
Kravtsov et al. 1998). Whether MOND is consistent with 
the RCs of dwarf galaxies, is an issue that is currently still 
under debate (Lake 1989; Milgrom 1991; Sanchez-Salcedo 
& Hidalgo-Gamez 1999). 

Recent studies have shown that more massive LSB disk 
galaxies reveal rotation curves that are similar to those of 
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dwarf galaxies (e.g., de Blok, McGaugh & van der Hulst 
1996; van Zee et al. 1997; Pickering et al. 1997), suggesting 
that the problem with the density profiles of dark matter 
halos is not limited to low-mass systems. Indeed, MB98a 
have argued that the observed rotation curve of LSB galax- 
ies are inconsistent with the NFW profile. MOND, how- 
ever, has been shown to yield remarkably good fits to the 
RCs of LSB galaxies (de Blok & McGaugh 1998; Sanders 
& Verheijen 1998). It is important to note, however, that 
the majority of the data on the RCs of the more massive 
LSB galaxies is severely affected by beam smearing (see 
e.g., Swaters 1999). Van den Bosch, Robertson & Dalcan- 
ton (1999) have shown that once these effects are taken 
into account, the RCs of the more massive LSB galaxies 
are in excellent agreement with dark matter halos with 
steep central cusps. 

5.1. A detailed comparison 

Given the unique opportunity to discriminate between 
DM and MOND based on the RC shapes of disk galax- 
ies, we now compare the RCs predicted by our DM and 
MOND models in some detail. Using our models which 
best fit the X-band TF relation (L5 and Ml), we have 
selected model galaxies that are photometrically similar 
(same scale length, same absolute magnitude, and same 
central surface brightness), and computed their RCs. In 
Figure ^| we compare four sets of model galaxies that span 
four magnitudes in both surface brightness and mass (see 
Table || for the parameters of each of the eight galaxies). 

For the high surface brightness galaxies, (dl, ml, d3, 
and m3) the RCs for the MOND and the DM models are 
very similar. For the LSB galaxies, however, the galax- 
ies from the DM model predict RCs that are significantly 
steeper in the inner parts. This is the reason why dark 
halo fits to several LSB galaxies have been claimed to fail, 
and owes to the steep central cusp (r _1 for the NFW pro- 
file used in our DM models) of the dark halo. This also 
explains why T(i?) is slightly overestimated for accelera- 
tions at around 3 x 10 9 <; a <, 10 10 ms~ 2 (see Figure 
as compared to the observations. Note, however, that the 
data on which Figure || is based suffers strongly from the 
beam-smearing effects eluded to above, and it remains to 
be seen whether it is actually inconsistent with the predic- 
tions from our DM model (see van den Bosch et al. 1999). 

The steepness of the central part of the RCs can be 
quantified by the parameter -R3/4, which is defined as the 
radius at which the circular velocity equals 75 percent of 
Vfiat- MB98a obtained R3/4 for a number of disk galax- 
ies with reasonably well resolved HI rotation curves, and 
found R^/^/Rd to strongly increase with decreasing lu- 
minosity and surface brightness. They compared this to 
predictions based on the models of Dalcanton et al. (1997; 
hereafter DSS97), which are similar to the models pre- 
sented here. According to MB98a, the models of DSS97 
predict an opposite behavior from that observed, with 
i? 3 / 4 /i?d decreasing with decreasing magnitude, and they 
use this as strong evidence against DM models. 

To further test this, we compute R3/4/ ' Rd for the model 
galaxies in each of our three models. Results are plotted 
as function of Mb and /io,B in Figure [h]. The first charac- 
teristic to notice is that none of the models reveals a sig- 
nificant increase of R 3 ^/Rd with decreasing magnitude. 
The MOND models reveal a modest increase of R^/^/Rd 



with decreasing surface brightness, but they reach a max- 
imum at -R3/4 ~ 1.25Rd- The galaxies of model L5 never 
reach values of R^/^/Rd in excess of ~ 0.8. Both mod- 
els are therewith in clear contradiction with the observed 
RCs presented by MB98a, which reach values as high as 
R^/i/Rd ~ 3. The most likely explanation for this dis- 
crepancy is, once again, beam smearing, which tends to 
over-predict Rz/^/Rd (see e.g., Blais-Ouellette, Carignan 
& Amram 1998). Data of higher spatial resolution is re- 
quired to test this further. However, given the modest 
difference between the predictions from MOND and the 
DM models, it seems unlikely that R^/^/Rd will allow to 
discriminate between the two scenarios 

Not only are the results presented here inconsistent with 
data, they also disagree with the curves plotted by MB98a 
which they claim to represent the model predictions from 
DSS97. However, we do not understand how MB98a have 
obtained these predictions. Part of the answer may lie in 
the fact that MB98a quote that they used models with 

1 /3 

L cx S ' "as predicted by DSS97". However, the mod- 
els of DSS97 predict £0 L 1 ^ 3 . Nevertheless, even with 
this erroneous interpretation of DSS97's models, we are 
unable to reproduce the predictions of MB98a. The re- 
sults presented here suggest that neither MOND nor the 
DM hypothesis predict a strong dependence of R^/i/Rd 
on either luminosity or surface brightness. Together with 
the fact that the MOND and DM models yield results 
that differ only modestly, we conclude that they both are 
equally (in)consistent with the data. The observed values 
of Rz/4/Rd can not be used to simply rule against the DM 
hypothesis only. 

6. CONCLUSIONS 

We have presented detailed models for the formation of 
disk galaxies in both a CDM and a MOND universe. In 
the case of DM, the structure of the disk is governed by 
the mass and angular momentum of the proto-galaxy. For 
MOND, however, the distribution of angular momenta of 
proto-galaxies is not known, and we have instead assigned 
scale lengths to the disks that are in agreement with ob- 
servations. In addition to these recipes that determine the 
structure and dynamics of the disks, we include recipes 
that describe how (part of) the gas is transformed into 
stars over the lifetime of the galaxy. We take account of a 
stability related threshold density for star formation and 
feedback from supernovae. The models have been tuned to 
fit the observed iC-band TF relation, and compared to nu- 
merous observations. Both models do remarkably well in 
reproducing a wide variety of observations of disk galaxies 
that span several orders of magnitude in both luminosity 
and surface brightness. 

To more easily compare the different models a summary 
of the results is given in Table ||, where we indicate, for 
each of the observational constraints discussed in this pa- 
per, whether the models are consistent with the data or 
not. In the case of DM (model L5), SN feedback is re- 
quired to yield a TF slope as steep as observed. Although 
the MOND model without feedback (Ml) is in excellent 
agreement with the empirical TF relation, it does not re- 
produce the observed deficit of HSB dwarf galaxies with 
small mass-discrepancies. This can be remedied by intro- 
ducing SN feedback (model M2), but at the cost of a TF 
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relation which is too steep and reveals an amount of scat- 
ter that is only marginally consistent with the data. This 
is a serious problem for MOND; since pure dynamics al- 
ready predicts a TF relation as steep as observed, it leaves 
virtually no room for any other galaxy characteristics to 
vary systematically with mass. 

The DM model is consistent with all observational con- 
straints against which we have tested it. Consequently, we 
strongly disagree with the picture that emerges from the 
literature (for example, MB98a, MB98b, and McGaugh 
1998), that the DM hypothesis suffers from numerous se- 
rious fine-tuning problems, which do not seem to have a 
clear-cut solution, whereas MOND is free from such prob- 
lems and capable of fitting virtually everything. We have 
shown here that once the DM model is tuned to fit the 
slope of the TF relation, it automatically passes all the 
tests devised by MB98a and McGaugh (1998) to argue 
against it. In particular, our DM model reproduces 

• A close correlation between global mass-to-light ra- 
tio and surface brightness, such that HSB and LSB 
galaxies follow the same TF relation without a sys- 
tematic offset. 

• The presence of a characteristic acceleration as ob- 
served. 

• A close relation between the characteristic acceler- 
ation, £, and central surface brightness of the form 

This is a remarkable result: there is no obvious reason 
why the DM model would reveal a characteristic acceler- 
ation, unlike in the case of MOND, where it is integral to 
the theory. Furthermore, it is encouraging that the same 
feedback parameters that yield the correct TF slope, result 
in an amount of scatter around the £-£o relation that is 
in excellent agreement with observations, and in addition 
explains the observed absence of HSB dwarfs with small 
mass-discrepancies. 

Of the list of observational facts about disk galaxies pre- 
sented in jy, the upper limit to the observed surface bright- 
ness of disk galaxies (item 5) has not been addressed by 
our models. Several studies have shown, however, that in 
the DM picture the presence of a maximum central surface 
brightness of disk galaxies is related to stability arguments 
(DSS97; Mo et al. 1998; Scorza & van den Bosch 1998). In 
§ |3.3| we have shown that the same stability argument also 
yields an upper limit on the central surface brightness of 
disks under the hypothesis of MOND (cf. Milgrom 1989). 
Henceforth, both the DM and the MOND models are con- 
sistent with the observational constraint of item 5, as long 
as disk stability is taken into account. 

The main problem for CDM is related to rotation 
curve shapes. We have compared predicted rotation curve 
shapes within the context of both MOND and CDM. Only 
in the case of LSB systems do the two scenarios yield 
RCs that are significantly different. Several studies have 
pointed out that dark halos with a steep cusp are incon- 
sistent with the central rotation curves of LSB and dwarf 
galaxies. This problem is also evident from the fact that 
at accelerations of ~ 1CP 10 hjo ms~ 2 the DM model pre- 
dicts mass-to-light ratios that are slightly too high. How- 
ever, it is important to realize that most data on the more 



massive LSB disks is severely affected by beam-smearing. 
This tends to underestimate the central gradients of ro- 
tation velocities, especially in galaxies that have a cen- 
tral hole in their HI distribution. When beam smearing 
is taken into account, the HI rotation curves of massive 
LSB galaxies are consistent with dark halos that follow a 
NFW density profile (van den Bosch et al. 1999). Further 
indications that beam smearing plays an important role 
comes from a comparison of the ratio R3/4/ ' Rd- Both the 
DM and the MOND models predict ratios that are lower 
than observed, which is most likely due to beam smearing. 
Therefore, in order to discriminate between MOND and 
DM, high resolution rotation curves of LSB disk galax- 
ies are required. Such data is currently only available for 
dwarf galaxies, which, because of their relative proximity, 
have been observed with high spatial resolution. It has 
been demonstrated convincingly that the rotation curves 
of these low-mass systems are inconsistent with centrally 
cusped dark matter halos. It is currently still under debate 
whether MOND can fit these rotation curves. If it can, this 
is where MOND has a clear advantage over CDM, unless 
(close to) constant density cores can be produced in dark 
halos that form in a CDM Universe (see e.g., Navarro, Eke 
& Frenk 1996, Kravtsov et al. 1998, and Bullock et al. 1999 
for possible solutions). 

One might argue that within the DM scenario one can 
fit basically anything as long as there are a sufficient num- 
ber of free parameters. In that respect it is important to 
realize that model L5 has only a very limited number of 
truly free parameters. Where possible, we have used pa- 
rameters that have either empirically determined values, 
or that are otherwise constrained: the dark halo proper- 
ties are taken from high resolution numerical simulations 
combined with the Press-Schechter formalism, T* K is con- 
strained by stellar population models, the star formation 
recipe uses values for Q and the Schmidt law that have 
been determined empirically, and a cr ;t is constrained by 
numerical simulations. This leaves only £g N and v as real 
free parameters. Tuning these two parameters to obtain 
a TF relation with the observed slope of b = —10.5, the 
model predicts gas mass fractions, characteristic accelera- 
tions, an To - So "conspiracy" , and global mass-to-light 
ratios which are in excellent agreement with observations, 
without additional tuning of the parameters. Further- 
more, we note that while MOND is often presented as 
being nearly free of fine-tuning, in order to match the sys- 
tematic properties of disk galaxies it is necessary to adjust 
the MOND feedback parameters to the same degree as 
required for DM. 

Probably the most amazing aspect of the models pre- 
sented here, is that the DM and MOND models are so 
very similar. However, both MOND and DM were con- 
structed to fit the rotation curves of disk galaxies. The 
fact that both theories correctly predict many other prop- 
erties, which are themselves closely related to the inter- 
nal dynamics (i.e., TF relation, gas mass fractions that 
are set by stability related threshold densities, SN feed- 
back whose efficiency depends on the escape velocity, etc), 
should therefore not be seen as too remarkable. It has 
often been argued that even if MOND turns out to not 
be correct, one should provide an explanation as to why 
it fits the properties of galaxies so well. The demonstra- 
tions in this paper suggest that any theory which yields 
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stable disks and fits their rotation curves, would probably 
perform as well as any of the models presented here. 
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APPENDIX 

A. THE B-K COLOR MAGNITUDE RELATION FOR SPIRAL GALAXIES 

In order to compare our models, which yield iVband magnitudes and surface brightnesses, to £?-band data, we require 
an estimate of B — K colors for disk galaxies. 

To determine an empirical B — K color magnitude relation we compiled apparent magnitudes in B and K as well as 
distances for a total of 139 spiral galaxies of type Sb or later from the following sources: forty-two Ursa-Major spirals 
from the sample of Verheijen (1997), sixty-one galaxies from the sample of de Jong & van der Kruit (1994), and thirty-six 
galaxies from the sample of edge-on spirals of Dalcanton et al. (1999). Absolute magnitudes were computed assuming 
H = 70km s~ 1 Mpc~ 1 and using the distances as quoted by these authors. Finally, we used the reddening maps of Schlegel, 
Finkbeiner & Davis (1998), with the R v = 3.1 extinction law of Cardelli, Clayton & Mathis (1989) and O'Donnell (1994), 
to correct for external extinction. 

The resulting B — K color magnitude is shown in the upper left panel of Figure [ll|. The solid line is the best linear fit 
with 

B - K = -0.66 - 0.18(M K - 5\ogh 70 ). (Al) 

The upper right panel plots the predicted -B-band magnitude, Mg rcd = 0.82M£- bs — 0.66 versus the observed value (for 

^70 = !)■ The data covers over eight magnitudes in Mb and reveals that M^ ed is a fairly good estimate of the observed 
B-band magnitude (with a standard deviation of <jm = 0.37 mag). 

The lower left panel of Figure [ll] plots the B — K color of the central surface brightness of the disks, versus /io.-ff for the 
data of Verheijen (1997) and de Jong & van der Kruit (1994). Unfortunately, at this time no central surface brightnesses 
are available for the data of Dalcanton et al. (1999). Note that B — K and ^o,s — Mo,i<" are generally not the same, since 
the disk scale lengths derived from the B and if -band data can differ considerably (see e.g., de Jong 1996b; Verheijen 
1997). The lower right panel of Figure [n] plots the predicted central surface brightness in the i?-band, ^'g, versus the 

observed value. Here is calculated from /J,Q h ^ and the color-magnitude relation of equation ( ]Al[ ). Note that the 

agreement is reasonable, but that the scatter is, with er M = 0.47 mag, slightly higher than for the absolute magnitudes. 

B. THE CHARACTERISTIC ACCELERATION PARAMETER £ 



In §4.4 we show that both the DM and MOND models yield a narrow correlation between the central surface brightness 
of the disk and the parameter £ defined by equation (fl3|). Here we investigate the reason for this narrow relation. 

In the case of MOND, VgL — G ao M ga i, and, if we assume that the galaxy is a pure exponential disk, it is straightforward 



to show that 



log (£mond) = ~0*o,b - 27.052) + \ log ( ^^T ga i ) , (B! i 

5 2 \ ao 



(see MB98b) , with T ga i the total mass-to- light ratio of the galaxy. 

In order to obtain a similar expression under the hypothesis of DM, we make the simplifying assumption that dark 
halos are represented by isothermal spheres, and that Vfl at = V200- Upon writing M4 — e g f /b ar M200 ( c f- equation Q), 
and since 

M *» = iJk' (B2) 

we obtain 



20 irGH V 200 Sp 

4dm - \ 5 J ■ ( a6 ) 

V a £ gf /bar 

If we now set Hq — 70 kms^ 1 and /bar = 0.085 (as for model L5, see Paper I), we can write 

log (4dm) = log (4mond) - 0.11 + ilog ( 250 ^° s -i ) - ^log(egf). (B4) 

Under the DM hypothesis one thus expects a similar relation between 4 and /io as under MOND, apart from an offset 

which depends on V200 an d the galaxy formation efficiency e g f . Since disks with the same central surface brightness can be 

embedded in halos with different rotation velocities (cf. Figure^), one expects a reasonable amount of scatter; a variation 

in V200 of a factor five, yields a scatter of 0.35 in log(£). However, this can be significantly reduced when SN feedback 

is important; the feedback efficiency will be larger, leading to smaller values of e g f in systems with lower values of V200, 

thus reducing the amount of scatter. Indeed, the scatter revealed by model L5 (see upper panel of Figure 0) is somewhat 

larger than in the case of MOND, but significantly smaller than 0.35 and therewith consistent with the data. 

1/2 

We thus conclude that the observed relation 4 oc S does not discriminate between the DM and the MOND hypotheses. 
However, the small amount of scatter observed may be considered evidence for efficient SN feedback, at least within the 
context of DM. 
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Fig. 1. — A comparison of the observed X-band TF relation with our three models. The upper left panel plots, with open circles, the twenty- 
two Ursa-Major spirals from the "unperturbed sample" of Verheijen (1997). The thick solid line is the best fitting TF relation (equation fet). 
This line is reproduced in each of the panels. The upper right panel plots the results for our best-fitting DM model (model L5), whereas 
the lower two panels correspond to the two MOND models discussed in the text. The value of the slope, b, of the best-fitting TF relation is 
indicated in each panel. 
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Fig. 2. — The HI mass-to-light ratios, Mm/Lg, as function of total magnitude. The upper left panel plots the data for a sample of one 
hundred disk galaxies of type Sb or later compiled by McGaugh & de Blok (1997), and clearly reveal an increase of gas mass fractions with 
decreasing luminosity. The remaining three panels plot the results for each of the three models, as labeled. The two thin lines have no physical 
meaning, but are plotted to facilitate a comparison between models and data. Luminosities and magnitudes for the model galaxies have been 
converted from the K-band to the £?-band using the B — K color magnitude relation presented in Appendix A. Given the errorbars associated 
with the data, and with the color conversion, all three models arc consistent with the data, and nicely reproduce the increase of Mm/Lg 
with decreasing luminosity. 



van den Bosch & Dalcanton 



15 



J 



m 



x 

130 
O 



© 



m 



3C 

tafl 
O 



1 - 







1 - 



- 




~ -1 - 



24 22 20 24 22 20 

Mo.b (mag/arcsec 2 ) yu 0B (mag/arcsec 8 ) 

Fig. 3. — Same as Figure^ except that now Mm/ Lb is plotted versus the central surface brightness of the disk, (j,o,b- Again, the thin 
lines have no physical meaning but are plotted to facilitate a comparison. Clearly, the three models are all consistent with the data and with 
each other. 
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Fig. 4. — The characteristic global mass-to-light ratio To (equation [113] ) as function of absolute B-band magnitude. The data (upper left 
panel) are taken from Verheijen (199 7; op en circles), McGaugh & de Blok (1998a; solid circles), and van Zee et al. (1997; solid squares). 
Details on the data can be found in j fe-U The data clearly reveal an absence of faint galaxies with low values of To, i.e., low luminosity 
systems always have a large mass discrepancy. This is reproduced nicely by model L5 (upper right panel), which is in reasonable agreement 
with the data (see the text for a more detailed discussion). Model Ml, the MOND model without SN feedback, however, is clearly inconsistent 
with the data, in that it predicts no deficit of low luminosity systems with low values of To. In model M2, we have included SN feedback, and 
tuned the corresponding parameters to reproduce the data. The thin solid line is plotted to facilitate a comparison, and is chosen to roughly 
outline the boundary of the data. 
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Fig. 5. — The central surface brightness of the disks^(/^o,s) as function of the rotation velocity at the flat part of the rotation curve, Va a t. 
The symbols for the data are the same as in Figure III The data reveal a clear absence of HSB galaxies with low rotation velocities. This 
is remarkably well reproduced by models L5 and M2 (both of which include SN feedback), whereas model Ml (no SN feedback) is in clear 
contradiction with the data. The thin solid line is chosen to roughly outline the boundary of the data, and has no further physical meaning. 
Note that the absence of observed galaxies in the upper left part of the diagram is not related to observational bias effects. Rather, it seems 
to indicate that SN feedback has played an important role in shaping low mass disk galaxies. 
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Fig. 6. — Same as Figure V&, except that we now plot the characteristic mass-to-light ratio To as function of the disk's B-band central 
surface brightness f^o,B- The oata reveal a narrow correlation, which is extremely well reproduced by each of the three models. The thin solid 
line corresponds to Tg = constant, and is plotted with arbitrary normalization. See text for a detailed discussion. 
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Fig. 8. — The ratio of total to luminous mass, T(i?) (equation ]14j ) . as function of radius R (panels on the left), orbital frequency uj (middle 
panels), and centripetal acceleration a (panels on the right). The upper panels plot the data for thirty galaxies from McGaugh (1998). 
These clearly reveal that the amount of scatter in T is minimized when plotted versus a, therewith indicating the presence of a characteristic 
acceleration. This character of the data is extremely well reproduced by the models. Whereas this is not surprising for the MOND models 
(Ml and M2), it is a remarkable success for the DM model. 
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Fig. 9. — Rotation curves for model galaxies plotted up to the radius at which the HI column density reaches 10 20 cm~ 2 . Thick solid 
lines correspond to the galaxies from model Ml (MOND), whereas the thick dashed lines are the rotation curves for galaxies from model L5 
(CDM). The thin dotted lines and long-dashed lines correspond to the contributions from the disk and the dark halo, respectively. Each panel 
compares the rotation curves of two galaxies that have virtually identical photometric properties (see Table 2), and which are indicated in 
each panel. The four sets of galaxies span four magnitudes in both luminosity and surface brightness. Note that the difference between the 
DM rotation curves and those for MOND increase with decreasing surface brightness. 
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defined as the radius where the RC reaches a velocity 75 percent of Vfl at , to the disk scale length as 
function of absolute B-band magnitude (upper panels) and central B-band surface brightness (lower panels) for each of the three models. 
Rotation curves in MOND are shallower (i.e., larger Rs/i/Rd) than in the case of DM, but only marginally so. Note that there is only a very 
marginal increase of R^/^/R^ with decreasing luminosity and surface brightness, in clear contradiction with the data presented in MB98a. 
In addition, both the DM and the MOND models, reveal a clear upper limit of R3/4 < 0.8Rd and -R3/4 < 1.25i?d respectively. This again is 
in clear contrast to the data, which shows galaxies with R 3 / 4 /R d as high as 3. The most likely explanation for this inconsistency seems that 
the data have not been properly corrected for beam-smearing (see discussion in text). 
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Fig. 11. — The upper left panel plots B — K colors versus absolute X-band magnitude for an ensemble of disk galaxies. Data is taken 
from Verheijen (1997; solid circles), de Jong & van der Kruit (1994; open circles), and Dalcanton et al. (1999; crosses). All galaxies have been 
converted to a common distance scale with Ho = 70kms — 1 Mpc -1 and have been corrected for external extinction in a consistent way, using 
the reddening maps of Schlegel et al. (1998). No correction for internal extinction has been applied to the data. The thick solid line in the 
upper left panel is the best-fit linear color magnitude relation (equation ]A1| ]). The upper right panel plots the predicted S-band magnitudes 
(based on the observed if-band magnitude and the best fitting B — K color magnitude relation) versus the observed values. The agreement 
is reasonable with a standard deviation around the line Mg Cd = M^ bs (thin line) of 0.37 mag. The lower two panels plot similar figures, but 
for the central surface brightnesses, rather than the absolute magnitudes (see Appendix A for further details). 
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Table 1 

Parameters of models discussed in the text. 



ID. 


Hypothesis 


Q 


1 K 


fc SN 


V 


b 


cm 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


L5 


CDM 


1.5 


0.40 


0.05 


-0.3 


-10.5 


0.20 


Ml 


MOND 


1.5 


0.53 


0.0 


0.0 


-10.6 


0.21 


M2 


MOND 


1.5 


0.53 


0.05 


-3.0 


-11.9 


0.33 



Note. — Column (1) lists the model ID, by which we refer to it in the text. Column (2) indicates whether the model 
is based on the DM or the MOND hypothesis. Columns (3) and (4) list the value of the Toomre parameter, Q, and 
the stellar if -band mass-to- light ratio in h^o M©/ L©. Columns (5) and (6) list the feedback parameters £g N and v, 
respectively. Finally, columns (7) and (8) list the slope b and scatter um (in mag) of the best fitting TF relation for model 
galaxies with -19 > M K - 51ogft, > -24 
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Table 2 

Parameters of model galaxies discussed in the text. 



ID. 


Model 


M K 


Rd 


Mo, if 


Vflat 


Ra/i/Rd 


A 


c 






mag 


kpc 


mag arcsec~ 2 


kms -1 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


dl 


L5 


-24.8 


3.8 


16.7 


237.5 


0.52 


0.040 


8.5 


ml 


Ml 


-24.8 


3.8 


16.7 


216.2 


0.68 






d2 


L5 


-24.8 


24.7 


20.8 


213.9 


0.41 


0.179 


8.3 


m2 


Ml 


-24.8 


25.0 


20.8 


233.8 


1.17 






d3 


L5 


-20.7 


1.4 


18.5 


92.1 


0.63 


0.030 


9.4 


m3 


Ml 


-20.7 


1.4 


18.6 


90.8 


0.98 






d4 


L5 


-20.7 


8.3 


22.4 


97.1 


0.46 


0.144 


9.8 


m4 


Ml 


-20.7 


8.3 


22.5 


114.9 


1.19 







Note. — Column (1) lists the model galaxy ID, by which we refer to it in the text. Column (2) indicates the model on 
which the galaxy is based. Columns (3)-(5) list the if-band magnitude, Mr, the disk scaledength, Rd, and the central 
surface brightness, Ho,k, of the model galaxies, respectively. Column (6) lists the rotation velocities measured at a column 
density of N[HI] = 10 20 cm~ 2 . The ratio of i£ 3 / 4 to Rd is listed in column (7). Finally, columns (8) and (9) list the spin 
parameter A, and the concentration parameter c of the NFW halo profile (before adiabatic contraction). Since we only use 
these parameters in the DM models, no values are listed for the MOND galaxies. All values listed correspond to h 70 = 1 
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Table 3 
Results. 



Observation 


item 


figure 


L5 


Ml 


M2 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


near-IR TF relation 


1 


1 


+ 


+ 




M m /L B 


4 


2,3 


+ 


+ 


+ 


T vs. M B 


3 


4 


+ 




+ 


T VS. Hq,B 


3 


6 


+ 


+ 


+ 


M0 : S VS. Vflat 




5 


+ 




+ 


£ VS. jUo,B 




7 


+ 


+ 


+ 


T(i?,w,a) 


G 


8 


+/- 


+ 


+ 


RC shapes 


2 


9,10 


+/-? 


+? 


+? 



Note. — Column (1) indicates the observational constraint. Columns (2) and (3) refer to the item number of the list 
of observational facts presented in and to the number of the figure in which a comparison of the models with the data 
is presented. Columns (4) to (6) indicate whether each of the three models discussed in this paper are consistent (+), 
inconsistent (— ), or only marginally consistent (+/— ) with the data. A question mark indicates that more data and/or 
work is needed before a definite answer can be given. 



